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ABSTRACT 

Hydrodynamic stability has been a longstanding issue for the cloud model of the broad line 
region in active galactic nuclei. We argue that the clouds may be gravitationally bound to the 
supermassive black hole. If true, stabilisation by thermal pressure alone becomes even more 
difficult. We further argue that if magnetic fields should be present in such clouds at a level 
that could affect the stability properties, they need to be strong enough to compete with the 
radiation pressure on the cloud. This would imply magnetic field values of a few Gauss for a 
sample of Active Galactic Nuclei we draw from the literature. 

We then investigate the effect of several magnetic configurations on cloud stability in axi- 
symmetric magnetohydrodynamic simulations. For a purely azimuthal magnetic field which 
provides the dominant pressure support, the cloud first gets compressed by the opposing ra- 
diative and gravitational forces. The pressure inside the cloud then increases, and it expands 
vertically. Kelvin-Helmholtz and column density instability lead to a filamentary fragmenta- 
tion of the cloud. This radiative dispersion continues until the cloud is shredded down to the 
resolution level. For a helical magnetic field configuration, a much more stable cloud core sur- 
vives with a stationary density histogram which takes the form of a power law. Our simulated 
clouds develop sub-Alfvenic internal motions on the level of a few hundred km/s. 

Key words: Galaxies: active, galaxies: nuclei, ISM: structure, hydrodynamics, radiative 
transfer 



1 INTRODUCTION 

Broad emission lines are produced in the immediate vicinity of op- 
tica lly active su per-massi ve black holes (SMBH, f or reviews see 
e.g. lOsterbrock 1988; Pe tersonlll997l : lNetzedl2008h . They may be 
used to infer the black hole mass in gala xies with such an Active 
Galactic Nucleus (AGN, e.g. lBentz et al.ll2009.) . and are a standard 
part of optically active AGN. The basic line emitting entity is usu- 
ally referred to as a cloud. The emission mechanism is photoionisa- 
tion by the central parts of the accretion disc. Photoionisation mod- 
els predict the clouds to have a typical temperature of order 10"* K, 
number densities of Wcl = lO''**' cm^-', sizes of Rc\ = 10^^^* cm 
and column densities of A'ci > 2 x 10 cm ' (e.g. lKwan & KrolikI 
1 19811 : iFerland&Elitzutll 19841 : iRees et al.lll989l) . Further important 
constraints come from reverberation mapping (e.g. jPeterson 1981, 
see below for a comparison of results of these two methods). The 
complete physics of these clouds is however highly complex and in- 
volves also pressure, radiative, centrifugal, gravitational, and mag- 
netic forces, probably on a very similar level. We are not aware of 
any attempt to include all these processes into a single model, but 
different authors have looked at some particular aspects of the prob- 
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lem: A general assumption for the cloud ensemble is often virial 
equilibrium. Since the gravitational potential is dominated by the 
SMBH, this would imply Kepler orbits. This treatment neglects the 
contribution of radiation pressure to the dynamics. The latter re- 
stricti on has been relaxed in a recent series of papers I Marconi et al 
,200^ l2009l : iNetzerl l2009l l2010l : iKrause. Burkert & Schartmann 
1201 ll (hereafter: paper I)). The general finding is that the radiation 
force may contribute significantly, and that in this case, the clouds 
on bound orbits could be significantly sub-Keplerian. For low angu- 
lar momentum (strong radiation pressure support), the orbits have 
to be highly excentric. 

The importanc e of radiation pressure has actually been recog- 
nised early on (e.g. 'Tarter & McKee"l973':'McKee & Tartedll975l: 
[Blumenthal & Mathews 1975, 1979; Mathews 1976). In fact, many 
of the earlier articles study the possibility that radiation pressure is 
dominant, also compared to gravity, and the clouds are unbound. In 
this case, it would however not be easy to understand why the black 
hole masses derived from reverberation mapping under the assump- 
tion of virial equilibrium may be brought into agreement with the 
correlation between black h ole mass and host galaxy properties 
with a uniform scaling factor dOnken et al.l2004l) . For optically thin 
clouds, one would expect almost no response to the variability of 
the ionising continuum. The strong variability of the emission lines. 
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lead by variations in the continuum is therefore evidence for tli e 
presence of optically thick clouds (e.g. ISnedden & Gaskelill2007h . 
However, the relative importance of gravity versus radiation pres- 
sure is proportional to the optical thickness (compare equation 
below). 

This opened up the possibility that the Broad Line Region (BLR) is 
gravitationally bound and possibly disc-like with a significant total 
angular momentum. In fact, several pieces of evidence that point in 
this direction have been found over recent years (compare paper I 
and references therein). One recent piece of evidence comes from 
spectropolarimetry, where the BLR is spatially resolved by an equa- 
torial scattering region, leading to differ ent polarisation an gles in 
the red a nd blue wings of emission l ines jSmith et al.ll2005h . Most 
recently, iKollatschnv & Zetzll ( l201lh have shown that the shape of 
the broad emission lines in many objects may be well fit with the 
assumption of a turbulent thick disc. 

Cloud stability and con finement is a long-standing issue 
jOsterbrock & Mathewsl 1 19861 for a review): The clouds should 
be in rough pressure equilibrium with th eir environment (e.g. 
iKrolik, McKee. & TarteJl98ll : lKrolikll988h . To reach the required 
pressure of p ~ 10^^ dyne cm^^, the inter-cloud medium needs ei- 
ther a high tem perature (> 4 x 10^ K, Krolik 1988) or a strong 
magnetic field jReej [l987h . Apart from the c onfinement issue , 
the clouds should be hydrodynamically unstable: iMathewsl il982l) 
has shown that while optically thin clouds may be close to uni- 
formly accelerated by the radiation pressure, there remain inter- 
nal radial pressure imbalances, which especially in the optically 
thick case, which is preferred by photoionisation models, lead to 
lateral expansion of the clouds. He termed the latter state pancake 
clouds. IMathewsl ( Il986l) then assessed the hydrodynamic stability 
of such pancake clouds with the result that the lateral edges of a 
pancake cloud are hydrodynamically unstable. Hence, the lateral 
flows persist and destroy the cloud on a short timescale, compa- 
rable to one cloud orbit. Also, the ram pressure by the inter-cloud 
gas, which is probably pressure supported and moving in a dif- 
ferent way than the clouds, may compress the clouds in the di- 
rection of relative motion and the incr eased pressure makes the 
clouds expand sideways iKrolikl f 19881) . now with regard to the 
direction of motion. These problems have le ad to the idea that 
clouds must reform or be re-injected steadily. iKrolikI l ll988l) de- 
velops the idea of clouds forming by the thermal instability from 
the hot inter-cloud gas. T his idea has however been rejected by 
iMathews & Doand jl990l) : It would require very large amplitude 
fluctuations of unusual type (high density, low temperature). Ad- 
ditionally, a compression by a factor of 10,000 would require an 
unusually small magnetic field, for the ma gnetic pressure not to 
stop the collapse. iMathews & Doan3 \\99& i further argue that the 
clouds would be destroyed by the radiative shear mechanism be- 
fore they could contribute to the emission line profiles. We have 
studied the radiative sheari ng mechanism for dusty clouds in 2.5D 
hydr odynamic simul atio ns jSchartmann et~^ 201 1). In agreement 
with lMathew3 ( Il986l) and lMathews & Doanj ( ll99QD . we find quick 
radiative shearing of the cloud. 

As mentioned above, the stability considerations were for 
clouds accelerating due to the dominant radiation pressure. The 
problem is even more severe for bound clouds where gravity is 
comparable to the radiation force: The radiation force acts on the 
illuminated surface of the cloud, while every part of the cloud is 
uniformly subject to gravity which must dominate by definition 
for bound clouds. Such clouds are therefore compressed, and must 
consequently be stabilised by some internal pressure. Yet, radia- 
tion pressure and gravity act in one dimension, whereas the thermal 



pressure is isotropic. It is therefore not possible to stabilise illumi- 
nated bound clouds by thermal pressure alone. 

In this context, we investigate the stability of magnetised and 
gravitationally bound BLR clouds - initially close to an equilib- 
rium orbit as calculated in paper I. Magnetic fields have so far 
rarely been considered in BLR clouds. We therefore first review 
literature data for indications on the relative magnitudes of grav- 
itational, magnetic and radiative forces, (section [2}. We then fo- 
cus on the effect of the magnetic field, and simulate the evolu- 
tion of irradiated magnetically dominated bound clouds. We study 
the two-dimensional, axisymmetric (2.5D), magnetohydrodynamic 
(MHD) evolution of isolated, initially optically thick clouds, in- 
cluding gravity, rotation and radial radiation pressure via a simple 
equilibrium photoionisation ansatz. We neglect self-gravity, viscos- 
ity and any radiation source other than the central accretion disc. 
Because of the small size of the clouds, our computational domain 
is small compared to the full size of the BLR. We describe techni- 
cal details in section[3] the setup details in section|4l our results in 
section |5] and discuss our findings in section (6] We summarise our 
results in section [7] 



2 MAGNETIC FIELDS, GRAVITY AND RADIATION 
PRESSURE IN BLR CLOUDS - INDICATIONS FROM 
THE LITERATURE 

The typical pressure level imposed on the clouds by radiative, and 
gravitational forces may be calculated from reverberation mapping 
(RM) data as follows: For this analysis we convert all forces to 
pressures, dividing by the surface area of the cloud, assumed to be 
spherical in this order of magnitude analysis. The inwards pressure 
due to gravity on the cloud is given by: 



PC ■- 



1 



(1) 



where Mgn is the black hole mass, m^xi the mass of the cloud, r the 
distance of the cloud from the black hole, and R^\i the radius of the 
cloud. In current reverberation mapping studies (e.g. iBentz et aEI 
[2009), the black hole mass is determined as a function of the mea- 
sured quantities of the centroid time lag Tcent and the line dispersion 
of the rest frame rms-spectra, dwa^ - Using these variables, the grav- 
itational pressure on the clouds may be expressed as: 



PC = ^f— PcldCline- 

^ ^cent 



(2) 



Here, p^ia is the density of the clouds, and / = 5.5 is the correction 
factor which ensures that RM-based black hole masses agree with 
the MBH-f* relation (Onken et al. 2004). We have collected Tcent 
and Oiine from a number of RM-studies where they were explicitly 
stated (TablefTJ. Where available, we have included different mea- 
surements of the same source at different epochs. From this data 
we have calculated pQ using /?cldPcld = 0.23 g cm^^ (correspond- 
ing to a hydrogen column density of Nh = 10^"' cm^^), which is 
plotted against black hole mass in Figure [T] (top). The values for 
i^cld and Peid are chosen to be in the range allowed from photoion- 
isation models (compare section [T]!. There seems to be no obvious 
correlation. The mean of the gravitational pressure in the clouds 
is 3.63 dyne/cm^. Because the distribution seems to be more uni- 
form in log-space, we use in the following the median value, which 
is 1.75 dyne/cm^, as a characteristic number. 

The ratio between radiation pressure and gravitational pres- 
sure may be expressed as: 
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l67tcG i?cldPcld Mbh ' 
where Ljo,, is that part of the AGN hght thatis absorbed by the 
cloud. It is given by (e.g. iMarconi et al. lEOOi) Lion = flfccLs 100, 
where a = 0.6 is the ionising fraction of the bolometric luminos- 
ity, b = 9 is the bolometric correction and L5100 represents XL^ at 
5100 A. 

The factor c accounts for geometrical and relativistic ef- 
fects: ISun & Malkarj ( Il989l) have shown that for frequencies up 
to about V = 10^^ Hz, the angular radiation pattern of thin ac- 
cretion discs around fast rotating black holes follow the classi- 
cal cos(9)-law (Lambert's law). For v > lO"" Hz, the emission 
is close to isotropic. For Schwarzschild black holes, the tran- 
sition occurs at somewhat higher frequencies. Luminous AGN 
emit most of their energy in the innermost regions of the accre- 
tion disc, in the blue - ultraviolet part of the spectrum (so-called 
"big blue bump"), with a peak arou nd Vp = 2 x lO'^ Hz (e.g. 
iHopkins. Richards. & Hernauistl20()7l) . If Lambert's law would ap- 
ply and if the BLR would extend 20 degrees above and below the 
equatorial plane, which would mean that about 10 per cent of the 
AGN's radiation would be intercepted, c would average to 0.3. For 
more rapidly spinning black holes or accretion discs emitting at 
higher frequencies, c might become as large as unity. Also, if the 
BLR would not be located in the equatorial plane of the accretion 
disc, c would also tend towards unity. Here, we adopt c = 1 as pos- 
sible, reasonable value. 

For these ass umptions, we have ca lculated the pressure ratio 
for the sample of IPeterson et alj ( |2004|) (their complete Table 8). 
We have used this sample, because black hole masses determined 
in the same way as for the sample above as well as continuum lu- 
minosities were available. We show a plot of the derived pressure 
ratios against black hole mass in Figure [T] (bottom). There is again 
no obvious correlation. The mean ratio of radiation to gravitational 
pressure is 0.74. Here, the distribution is also more uniform in log- 
space, and therefore again we use the median of the pressure ratios 
as a characteristic number, which is 0.48. These numbers depend 
of course crucially on the assumed column density. Photoionisa- 
tion calculations constr ain the column density to A' > 2 x 10^^ cm^^ 
dKwan & Krolik|[l98 lV We have used values for and p^id cor- 
responding to a (hydrogen) column of N = lO^'^cm^^, consistent 
with the requirement from reverberation mapping studies that the 
clouds need to be optically thick (compare section [T] above). Al- 
lowing fo r the smallest column density consistent with the calcula- 
tions of Kwan & Kroliklh98lh would shift the median of the pres- 
sure ratio to 2.39, w i th onl y 6 o f 35 values st i ll belo w unity. How- 
ever, |M^oljno|eral] j20ICl) and lRisaliti et alj ( l201lh have possibly 
observed BLR clouds in X-ray absorption and give values of a few 
times lO^-' for the column of their clouds. With such a column den- 
sity, the typical pressure ratio would clearly drop below unity. We 
have included a contribution of 80 per cent due to the uncertainties 
in the column density for the calculation of the error bars in Fig- 
ure[T] Another uncertainty is due to the accuracy of the black hole 
masses: As indication, the errors stated by Bentz et al. (2009, their 
Table 13) ar e 43 per cent. Add itionally, there may be a systematic 
error via / jOnken et al.ll2004), and possi bly corrections to radi- 
ation pressure ( Marconi et alj 20081 . 1 20091) . Because the pressure 
ratio values we find scatter around unity (Figure [TJ, it is thus not 
possible to firmly conclude if the BLRs are typically gravitation- 
ally bound or otherwise. It is however remarkable that the pressure 
ratio distribution is bounded by a value close to unity for reason- 
able assumptions. In any case, gravitationally bound cloud models 
are clearly consistent with the data. 
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Figure 1. Radial pressure on BLR clouds. Top: Inwards pressure due to 
gravity over black hole mass (data of Table [T). Bottom: ratio o f radiative 
to gravitational pressure (all 35 objects of Pete rson et alj|2004l) . Different 
samples have been used for each plot, according to the availability of the 
data. Error bars are propagated from the original publications (compare Ta- 
ble[T). Additionally, an error of 80 per cent on the assumed column density 
has been included. See text for details. 



The typical radiation pressure on BLR clouds is given by the 
characteristic gravitational pressure times the typical pressure ratio, 
and using the fiducial numbers derived above: 



PR = 0.84 dyne cm 



(4) 



This number exceeds the upper bound for the thermal pressure de- 
rived from standard photoionisation models (compare above) by a 
factor of a few. It is therefore non-negligible. Even, if one takes the 
lower bounds of the distributions of pc and Pr/pg, the resulting 
radiation pressure is well within the range of inferred values for the 
thermal pressure. 

The fact that the reverberation mapping results imply small 
BLR radii and therefore hi gh radiat i ve flu x es has already been 
noted by several authors jPetersorJ Il988l : iFerland et alj Il993 : 
iLeighlv & Casebeeill2007l) . It would seem to imply that the pho- 
toionisation parameter is high. Adopting S = 2.3pr/pt as pho- 
toionisation parameter ('Krolik et alj|l98lj) . we arrive at a charac- 
teristic value of E = 14, taking the upper bound of the thermal 
pressure from the photoionisation mo dels, and ten times higher for 
the central value. lKroliketalJ ( ll98lh point out two arguments to 
constrain E: First, pressure balance with a surrounding hot phase 
implies 1/3 < E < 10. Second, fitting the line ratios, requires 
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Table 1. Reverberation mapping data 



Object 






AfJi 




(days) 


(days) 


(km/s) 


(km/s) 




Bentz et al. 


. (2009) 






Mrkl42 


2.74 


0.75 


859 


102 


SBS1116+583A 


2.31 


0.56 


1528 


184 


ArplSl 


3.99 


0.59 


1252 


46 


Mrkl310 


3.66 


0.60 


755 


138 


Mrk202 


3.05 


1.43 


659 


65 


NGC4253 


6.16 


1.43 


516 


218 


NGC4748 


5.55 


1.92 


657 


91 




4.18 


1.08 


4270 


292 


NGC6814 


6.64 


0.89 


1610 


108 




Dennev et al. (2009) 






NGC4051 


1.87 


0.52 


927 


64 




Grier et al. 


(2008) 






PG2 130+099 


22.9 


4.7 


1246 


222 




Dennev et al. (2006) 








3.73 


0.75 


1 S^»l 
1 JOI 






Bentz et al. 


. (2007) 






NGC5548 


6.3 


2.5 


2939 


768 




Bentz et al. 


. (2006) 






NGC5548 


6.6 


1.0 


2680 


64 




Peterson et al. (2004) 






Mrk335 


16.8 


4.4 


917 


52 


Mrk335 


12.5 


6.1 


948 


113 


PG0026+129 


98.1 


26.9 


1961 


135 


PG0026+129 


111.0 


26.2 


1773 


285 


PG0052+251 


163.7 


48.4 


1913 


85 


PG0052+251 


89.8 


24.3 


1783 


86 


PG0052+251 


81.6 


17.7 


2230 


502 


Fairall9 


17.4 


3.8 


3787 


197 


Fairall9 


29.6 


13.7 


3201 


285 


Fairall9 


11.9 


5.7 


4120 


308 



Centroid time lag. 

Error on centroid time lag. 

Line dispersion, rest frame rms-spectra. 

Error on line dispersion. 



0.3 < E < 2 (lKwan&Krolikll98lh . Additionally, very high ionisa- 
tion parameters lead of course to highly ionised gas, which would 
be unable to produce the observed emission lines. Yet, the uncer- 
tainties in the column density, geometrical and relativistic correc- 
tion factors and black hole masses are considerable. Thus, while 
the ratio of radiation pressure to thermal pressure comes out a little 
high, which might point e.g. to a lower value of our factor c (com- 
pare above) or cloud densities which are generally rather around 
lO" cm^-^, the data still seems to be consistent. 

In summary, the data suggests that gravitational and radiation 
pressure are both important and at least comparable to the thermal 
pressure. Gravitationally bound clouds are are consistent with all 
the available data. 

We have argued above that thermal pressure, due to its 
isotropic nature may not stabilise the cloud against the uni- 
directional opposing forces due to gravity and radiation. Magnetic 
forces are not isotropic. Their direction depends on the geometry 



of the magnetic field. It is therefore natural to ask if a magnetic 
field configuration may be found that is able to support an illumi- 
nated bound cloud. If the magnetic field should be at all capable to 
contribute, its energy density should be at least comparable to the 
radiation pressure. Both, magnetic pressure and magnetic tension 
are of the same order of magnitude as the magnetic energy density 
for many, not highly symmetric configurations. Accordingly, the 
magnetic field strength required for the magnetic energy density to 
match our reference radiation pressure (equation {4}) is 5 G. Taking 
into account the accuracy of the measurements, as well as the width 
of the distribution of pressures (FigureO, the true indicative value 
for any given BLR should be within a factor of ten of this number. 

Guided by these considerations, we make the following as- 
sumptions for our magnetohydrodynamic stability analysis: 

(i) The clouds are gravitationally bound. 

(ii) The magnitude of the radiation pressure is significant com- 
pared to gravity. 

(iii) The magnetic pressure is at least comparable to the radia- 
tion pressure. 

(iv) The thermal pressure of the cloud is much smaller than both 
the radiation pressure and the magnetic pressure. 

The last assumption is mainly for methodological reasons, in or- 
der to isolate the effects of the magnetic field. In the following we 
present magnetohydrodynamic simulations of clouds set up accord- 
ing to the standard cloud properties outlined in section [T] and the 
above assumptions, and investigate their stability numerically. 



3 NUMERICAL MODEL 

The basic code we use is the 3D MHD code NIRVANA 
jZiegler &~ Yorke 1997). For short, it conserves mass, momentum 
and internal energy in the advection step , interpolates the fluxes 
with van Leer's formula Ivan Leeilll977l , second order accurate), 
and uses the constrained transport method to keep the magnetic 
field divergence free. 

In order to include the effects of photoionisation and electron 
scattering, we augmented the MHD equations by an equation for 
the radial radiative transfer for the photon flux S: 



dS . (2 , ^ 

= -"HI ^ I ■;: + ^T"Z ) s, 



(5) 



where the first term on the right-hand side is due to photoionisation, 
and the second one describes the geometrical effect and electron 
scattering, with the Thomson cross section Ox- In equilibrium, the 
ionisation rate nei is given by the number of recombinations: 



"HI = arec(«Z) 



(6) 



where n is the total number density of ionised and neutral hy- 
drogen atoms. We assu me a recombination coeffic ient of cHec = 
5 X 10"'^ cm^ s"' (e.g. lDopita & Sutherlandll2003h . kept constant 
for simplicity. The ionisation fraction % is given by: 

Cphof^ 



y 



ttrec" 



X ^ - 



(7) 



(8) 



We assume a photoionisation cross section of C7phot = 6.3 x 
10^'**/(a + 3.5) cm^, averaged over a spectrum v^" with a = 2. 
This leads to a local radiative acceleration of: 
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«rad = («HI+p!C!'T5)/lV/pc, 



(9) 



Shell position 



where p denotes the mass density and c the speed of hght. We as- 
sume an average photon energy hv = 1 .007^ Vg, with the hydrogen 
ionisation potential hVo = 13.6 eV, for a steeply falling spectrum. 
The radiative heating is not taken into account explicitly. We set 
however a lower limit to the radiative cooling at 1000 K, which im- 
plies a heat term that balances cooling in such a way as to reach the 
minimum temperature. We do not attempt to model the tempera- 
tures within the clouds correctly. This is justified, because we study 
the evolution of magnetically dominated clouds, and even if pho- 
toionisation heating was taken into account, the thermal pressure 
would remain negligible in our setup. We set the minimum temper- 
ature below the one expected from photoionisation equilibrium, in 
order to investigate if shock heating can also heat the clouds to a 
temperature comparable to the photoionisation temperature. 

Equilibrium photoionisation is a good assumption, because 
the ionisation front would pass the clouds quickly (< hour), and the 
thermal pressure imbalance induced by the passing of the ionisation 
front would be insignificant compared to the magnetic pressure. 
The clouds we consider below are about halfway ionised. Thomson 
scattering dominates in the inter-cloud gas, hydrogen photoionisa- 
tion within the clouds, which also have a still significant contribu- 
tion by Thomson scattering. This model was chosen with the aim 
to be as simple as possible, while capturing the essential physics 
(we are interested in the radiative acceleration). We have varied the 
parameters within reasonable limits without a qualitative change of 
the results. We adopt an equilibrium solar metallicity cooling curve 
( [Sutherland & Dopita 19933, extrapolate d to lower temperatures as 
described in lKrause & Alexandeij i l2007h . Finally, we take into ac- 
count the gravitational acceleration g due to a central point source 
(the SMBH). 

These effects are implemented as source terms into the system 
of MHD equations: 



dp 
dt 
dpv 



dt 



+ V.(pv) 
-V - (pvv) 

+ V-(ev) 
'dt 





- Vp + — (B • V)B - — VB^ 

+P(«rad-5) 
V • V - <^ 

Vx (vxB), 



(10) 

(11) 
(12) 

(13) 



where p denotes the density, e internal energy density, v velocity, 
the radiative cooling and p = (7— l)e the pressure. 
The problem we are trying to address here has some par- 
ticular requirements: We have verified that angular momen- 
tum is conserved to high accuracy during our simulation runs. 
iKrause & AlexandeJ ^2007,) have shown that energy is conserved 
well in multi-phase setups up to density contrasts of about 10^. 
Also in the simulations presented here, we have checked that the 
total energy is reasonably well conserved. We are however limited 
in density contrast to about 10^. 



4 SIMULATION SETUP 

We simulate isolated clouds in axisymmetry. The initial geometry 
is toroidal with a circular shape in the meridional section. The basic 
simulation parameters are typical for BLRs and are summarised in 
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time / yr 
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Figure 2. Radial position of the centre of mass of a ID test cloud (green 
plus-symbols) over time compai'ed to the analytic prediction (solid black 
line) for an optically thick cloud. 



table|2] The cloud is always positioned at a distance of 9200 /xpcQ 
from the central SMBH. We resolve the cloud diameter of 6 /ipc 
with 30 cells in standard runs, and up to 90 cells in the highest res- 
olution ones. We set the initial cloud density to 4 x lO'^mp cm^-^, 
where mp is the proton mass. This results in a central hydrogen col- 
umn density of 5 x lO^^cm^^, except for run R47, which has twice 
this value. 

The setup follows essentially the scenario of paper I: We as- 
sume a hot atmosphere in approximate hydrostatic equilibrium, 
where the pressure follows a power law p{x) = po'^ ' with x = 
'"/''eld s = 2. Here, r^jj denotes the radial position of the cloud 
(distance from the SMBH) and po = GMbhPo(1 ^ T) /''cld^ is what 
the inter-cloud thermal pressure would be at that location, with 



10 OTp cm 



-3 



and the Eddington factor F = 0. 1 (/ in paper I). 



Here, we take radiation pressure due to Thomson scattering into ac- 
count. G is the gravitational constant and Mgn the black hole mass, 
which is generally assumed to be 10** solar masses. The density is 
given by p(r) = pox'^'. The magnetic field (compare below) in the 
inter-cloud medium is perpendicular to the computational domain 
and contributes initially about ten per cent to the total pressure, 
above hydrostatic equilibrium. Yet, because of the cooling the at- 
mosphere still develops an inwards flow over the simulation time. 

The pressure level might appear high, yielding an inter-cloud 
temperature of order 10^ K. This is however necessary since we 
would like to test here the effect of dynamically relevant magnetic 
fields primarily inside the clouds. In order to yield a smaller ther- 
mal inter-cloud pressure, and therefore a temperature of order the 
Compton temperature, which one might expect, one would have 
to assume a much stronger magnetic field also in the inter-cloud 
medium. We defer the investigation of this situation to future work. 

According to paper I, force equilibrium is reached for the 
cloud, if the rotational velocity in Kepler units is given by: 



1 



3F 



(14) 



where A' = 2R^i^Pc\tj/mp is a measure of the column density of the 
cloud. In paper I, we have found that for high column densities, the 
clouds should be in a minimum of the effective potential, and there- 
fore remain stable at their radial position. For low column densities. 
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Table 2. Simulation parameters 



Label 


comment 


„. a 














5r^ 


vl 


R34 


eq-cld 


9100 


9250 


32 


92 


3 


750 


460 


0.2 


85 


R35 


neq-cld 


9120 


9320 


32 


93 


3 


1000 


460 


0.2 


85 


R36 


R34-hires 


9100 


9250 


32 


92 


3 


1500 


920 


0.1 


85 


R47 


R34-big-cld 


9000 


9300 


64 


188 


6 


750 


460 


0.4 


94 


R60 


helical 


9150 


9250 


32 


92 


3 


1000 


920 


0.1 


85 


R61 


helical-hires 


9150 


9250 


32 


92 


3 


1500 


1380 


0.067 


85 


R62 


helical-lores 


9150 


9250 


32 


92 


3 


500 


460 


0.2 


85 


R63 


helical-big 


9100 


9300 


64 


188 


3 


2000 


1840 


0.1 


85 



" Inner boundary of the computational domain in r direction in /ipc. Note that the origin is not within the computational 
domain. 

^ Outer boundary of the computational domain in r direction in jlpc. 

Size of computational domain in 6 direction, symmetrical above and below the equator in units of 10^^ ;r. 

Corresponding approximate domain size in vertical direction in /ipc. 
^ Radius of the cloud in ^Jpc. 
f Number of cells in radial direction. 
^ Number of cells in meridional direction. 
^ Uniform resolution in /ipc/cell. 

' Initial cloud orbital velocity as a fraction of the Kepler velocity in per cent. 



the clouds should be either ejected, or be sent inwards on an excen- 
tric orbit. For our pressure profile, the critical rotation velocity is 
75 per cent of the Kepler velocit}0. All our clouds are set up on the 
stable branch of the equilibrium curve (compare Table[2ll. We found 
experimentally that we require a rotational velocity of 102 per cent 
of the value given in Equation[T4]to ensure initial force equilibrium. 
The difference is due to the small amount of Thomson scattering in 
the ambient medium on the way to the cloud, which reduces the 
amount of radiation that actually arrives at the cloud. 

For runs R34, R36 and R47, we use an azimuthal magnetic 
field, only. The plasma j3 = inp/B^ is initially set to 10 outside the 
cloud. Inside the cloud, we set the temperature to the minimum one 
(1000 K), and adjust the azimuthal magnetic field to yield pressure 
balance. This results in a magnetic field about three times stronger 
than in the inter-cloud medium, and a j8 = lO^'*. Run R35 uses 
a uniform field strength inside and outside the cloud again with 
j8 = 10 in the inter-cloud medium. This cloud is therefore initially 
underpressured. Runs R60, R61, R62 and R63 have an additional 
poloidal field (in the cloud, only). We have again j8 = 10 in the 
inter-cloud medium. Within the cloud, the azimuthal field declines 
with distance from the cloud centre as a Gaussian: 



=B^,.ic(l+3exp(-10i?V'^dd))> 



(15) 



where B^ ^,. is the azimuthal inter-cloud field, and capital R refers 
to the distance from the cloud centre. Rf.\^ is the cloud radius. The 
poloidal field is set up as closed field loop in a given meridional 
plane. The 3D structure would be helical. We initialise the radial 
component by 



„ «''cld d 
B, = ——^B^sme. 
rsmO do ^ 



(16) 



We have tried different values for a. The one we report here is 
a = 2 X 10^^. For this case, the initial peak values of the toroidal 



^ Note that eq. (6) of paper I contains a mistake. The correct formula for the 
critical rotation velocity should read: Vj? = ( 1 + 3/ 2s) ^ ' , with the pressure 
power law index .s = 2 in our case. 



and poloidal field components are 103 G and 15 G, respectively. 
The magnetic field strength is for all runs generally around 30 G 
(the larger initial field in the runs with helical magnetic field decays 
quickly to this value), to ensure magnetically dominated clouds, 
following the considerations in section |2] 

All boundary conditions are set to zero gradient for all vari- 
ables. We believe that the best test for our radiation pressure mod- 
ule is the approximate stability of the hydrostatic halo, mentioned 
above. We have also checked ID-cloud acceleration without gravity 
(Figure[2}. In this test, the cloud accelerates a bit more slowly than 
expected from simple radiative acceleration of an optically thick 
cloud. This is, because the cloud also has to work against the ram 
pressure of the ambient medium. 



5 RESULTS 

5.1 General dynamics 

A general overview of the evolution of our simulated magnetised 
clouds is shown in Figures [3] |4] and |5] Movies are provided with 
the online version. 



5.1.1 Pressure equilibrium 

We first describe the evolution for run R34 (Figure [3}, set up in 
initial total pressure equilibrium with an azimuthal magnetic field, 
only. The cloud is stationary at its radial position due to the matched 
radiative and centrifugal forces on the one hand, and the gravita- 
tional force on the other hand. This leads to a radial compression 
of the cloud (10.93 morb (= milli-orbits)). At the same time, the 
magnetic pressure in the cloud increases. As there is no opposing 
force vertically, the cloud expands upwards and downwards, much 
like an open tube of toothpaste squeezed i n the middle. This evo- 
lution matches the picture that iMathewsl l ll983) described as the 
formation of quasar pancakes. Towards the edges of the cloud, the 
column density is lower, and therefore, there is a net outward force. 
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Figure 3. Evolution of run R34. The top six images show the logarithm of the density for different snapshot times, indicated in milli orbits (morb) and months 
on the individual images. The bottom three images show for the final snapshot from left to right: Total pressure, magnetic field strength and magnetosonic 
Mach number for the velocity component in the meridional plane. In the cold cloud gas, the magnetosonic speed is almost the same as the Alfven speed. 



This can be clearly identified in the 10.93 morb density plot. How- 
ever, on the very edge, there is a certain amount of mixing. This 
will always happen to some extent, as the contact surface gener- 
ally has to be resolved by a few grid cells. Wherever the cloud 
gas mixes with the ambient gas, the cloud looses rotational sup- 
port and falls inwards. There is also some drag due to the inflow 
of the cooling ambient gas. A Rayleigh-Taylor instability due to 
the initial cloud acceleration is clearly seen in the middle of the 
cloud. The cloud continues to expand vertically by the tube of 
toothpaste mechanism, with lower column density regions being 
pushed outwards and mixing regions coming back inwards. By 



40.08 morb, the cloud has dispersed into three major fragments 
with lower density (m ixed) filaments falling inwards. The Kelvin- 
Helmholtz-timescale dChandrasekhaj 1 1 96 ih is comparable to the 
evolution time: 

A / 77 \ -1/2 / V \ -1 

Tkh = 1 month — — r- , (17) 

/xpc V 0.005/ V200kms-'/ ' 

where A is the wavelength, t] the density ratio and v the shear ve- 
locity. As expected, some filaments show the typical rolls of the 
Kelvin-Helmholtz instability. In the final snapshot (54.66 morb), 
all of the cloud material has mixed with the ambient gas to some 
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Figure 4. Same as Figure[3]for run R35. 
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Figure 5. Same as Figure|3]for run R63. Only part of the computational domain is sliown. 
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R34: density histograms 
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R35: density histograms 
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R63: density histograms 
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Figure 6. Density histograms for various runs, indicated in the individual titles and snapshot times (order: rder: solid black, dotted blue, dashed green, dot- 
dashed bright green, dashed red, solid orange), indicated in milli-orbits (morb) next to the respective lines. Mixing generally moves the upper density cutoff 
towards lower density (leftwards) for later times. The exception is run R63 for which the histogram converges at f « 9 morb and also continuously extends 
towards densities higher than the initial cloud density. 



degree, and falls leftwards, towards the SMBH. Mixing is a reso- 
lution effect in our simulations, but would also happen in reality 
at scales below our resolution. The mixing is further illustrated 
in the density histograms, which we show for similar snapshot 
times in Figure|6] Between t = 30.36 morb and the final snapshot 
t = 54.66 morb, the upper density cutoff drops by about a factor 
of five, indicating that the last cloudlet cores have been destroyed 
by mixing. Consequently, at this time, we also observe the cloud 
complex to start falling inwards as a whole. During the whole evo- 
lution time the entire simulated region stays well in pressure equi- 
librium (Figure |3] bottom left), the typical deviation being within 
ten per cent. Also, the magnetic field in the cloud complex re- 
mains close to its initial value of 31 Gauss (Figure[3l bottom mid- 
dle). Acceleration is mediated by the compression of the magnetic 
field. We expect therefore the velocities to be limited by the Alfven 
speed, Ca = Bj ~ 500 km/s in the cloud and ~ 6000 km/s in 
the ambient gas. Indeed, the observed velocities in the filaments 
are typically around 30 per cent of the Alfven speed (Figure [3] 
bottom right, we use the magnetosonic speed in the Figures to 
show that also in the inter-cloud gas, the velocities remain sub- 
magnetosonic). 



5.1.2 Underpressured cloud 

We check the dependency on the initial condition in run R35 (Fig- 
ure |4]l: Here the field strength is not higher inside the cloud, but set 
to the same value as in the ambient medium at the corresponding 
radius. The latter is still kept at j3 = 10. Also, everything else is as 
in run R34. Now, the cloud is first compressed isotropically due to 
the initial underpressure before the radial (with respect to the black 
hole) compression phase. It overshoots and oscillates in size a few 
times. The acceleration produces Rayleigh-Taylor-instabilities on 
the cloud surface (/ = 4.82 morb). Then the cloud gets compressed 
anisotropically in radial (with respect to the black hole) direction 
(/ = 7.23 morb), and suffers a similar filamentation process than 
run R34, above {t = 9.64 and t = 19.29 morb). This initial phase 
proceeds significantly faster than for R34. Between t = 19.29 morb 
and / = 37.37 morb, the remaining dense cloudlet cores are dis- 
persed (compare Figure |6j, and at f = 37.37 morb, the cloud is 
infalling towards the SMBH. The cloud, and respectively its frag- 
ments, is quickly compressed to the equilibrium magnetic field of 
3 1 Gauss, gets quite close to pressure equilibrium and acquires ve- 
locities of about 30 per cent of the magnetosonic speed (Figure |4] 
bottom plots). Thus, with the exception of the initial phase, the evo- 
lution is very similar to the one of run R34. 
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5.1.3 Cloud with helical field 

A significantly different result is obtained for a cloud with a he- 
lical magnetic field (R63, Figure |5). The cloud is initially not in 
force equilibrium, and first oscillates in size a few times. Simi- 
lar to run R35, this has caused some filaments by f = 3.02 morb. 
These filaments spread and disperse in much the same way as for 
the other runs. However, in contrast to the other runs, the cloud 
core remains remarkably stable, and is not significantly compressed 
by the radiative, centrifugal and gravitational forces. This situa- 
tion remains essentially unchanged until the end of the simulation 
{t = 43.55 morb), which we generally take to be the latest time be- 
fore a significant amount of cloud gas has left the grid. The final 
equilibrium is characterised by a total overpressure of a factor of 
1.3 in the cloud core, balanced by the magnetic tension force (Fig- 
ure [5] bottom left). In the outer filaments, the magnetic field is of 
a very similar magnitude to the other simulations. It is about ten 
per cent higher in the cloud core. The velocities are typically about 
60 per cent of the magnetosonic speed, only slightly higher than 
in the previous simulations and still sub-magnetosonic. Also the 
density histograms (Figure |6]( are distinctly different from the pre- 
vious runs: Whereas the histograms for the runs without a poloidal 
field component did not converge with time, with the upper den- 
sity cutoff moving to lower densities at the end of the simulation, 
the density histograms for run R63 stay essentially constant from 
at least f = 24 morb up to the end of the simulation at f = 44 morb, 
which corresponds to almost 1 /2 of the simulation time. The occu- 
pied volume is roughly given by 

AV 



dlogp 



(18) 



and also extends to densities exceeding the initial cloud density by 
a factor of ten. 

Close-ups of the cloud core at the beginning and at the end 
of the simulation are shown in Figure |7] The cloud has essentially 
the same size as in the beginning. The density is however no longer 
homogeneous, but displays a round filamentary pattern, following 
the cloud surface. The magnetic field geometry is essentially iden- 
tical to the one of the initial condition: a dominantly azimuthal core 
with a helix winding around. The peak value of the azimuthal field 
has however reduced by a factor of three, whereas the poloidal field 
strength remained almost unchanged. 

Interestingly, many of the filaments that the cloud sheds are 
essentially devoid of poloidal field. Due to the solenoidal condition, 
filaments shed by the cloud need to have a bidirectional magnetic 
field. We have observed this at 0. 1 G level near the cloud surface. 
But since they are often at the resolution limit, the poloidal field 
tends to cancel numerically, while the toroidal one is unaffected by 
this. Without poloidal field in the filaments, they disperse much like 
the filaments in the other simulations. Yet, sometimes filaments are 
shed which are thick enough to keep their poloidal field. Two such 
examples are present in Figure[7] Interestingly, they also adopt the 
field geometry of the parent cloud. This indicates that we may have 
found a magnetic field configuration for such clouds. It is clear that 
a poloidal field stabilises the cloud significantly. 



5.2 Resolution dependence 

We have re-simulated R34 at twice the spatial resolution (R36, Fig- 
ure [S] middle). The evolution is quite similar to run R34, with 
the exception that the dispersion of the cloud and cloudlets hap- 
pen faster. This may be seen from the more uniform appearance of 



the density distributions (Figure [8}, but also from the density his- 
tograms (Figure[6}, which are shown at the same timesteps: While 
the histograms at the corresponding times are generally very sim- 
ilar, the density cutoffs at r = 54.66 morb differ by about 0.2 dex. 
Thus, the simulation can be regarded as overall numerically con- 
verged, but there is increased turbulent mixing at small scales for 
increasing resolution. 

We have also studied the resolution dependence for run R63. 
For this study, we have used a smaller grid, because we have found 
that reflection of sound waves at the grid boundaries somewhat in- 
fluences the motion of the filaments. R60 is the re-simulation of 
R63 at the same resolution. R62 uses 50 per cent less cells on a 
side, R61 has 50 per cent more cells on a side. Logarithmic density 
plots of all three runs at / = 48.58 morb are shown in Figure|9] The 
filamentary systems are in general quite similar, though they have 
of course a finer structure at higher resolution. Yet, their extent and 
therefore kinematics is similar. The cloud core has drifted outward 
slightly for the two higher resolution runs. It is less dense and mov- 
ing inwards for the low resolution run. The resolution is in this case 
not enough to resolve the cloud core structure, and therefore, the 
clouds disperses faster. This is also evident from the density his- 
tograms (FigurefTOt. The low resolution run lacks cloud material at 
high densities, which the two higher resolution clouds do show. It is 
interesting that this is in a way opposite to the case without poloidal 
field. There, higher resolution had made the cloud disperse quicker. 
Here, high resolution is essential to preserve the cloud core at high 
density. Our results are thus robust under resolution changes. 



5.3 Scaling with cloud size 

The simulations are not scalable, as the source terms introduce a 
dimensional scale to the problem. In order to explore the scaling 
behaviour, we re-simulated run R34 with all scales doubled (R47, 
Figure [8] right). As expected, the evolution now takes longer, but 
not nearly twice as long as for run R34. Otherwise, we find a very 
similar behaviour: Again the cloud is first compressed radially, then 
expands vertically by the tube of toothpaste mechanism, disperses 
into a filamentary system and then moves inwards due to the mixing 
related loss of rotational support and drag by the cooling infalling 
ambient gas. 



5.4 Turbulent velocities 

We define the mass-weighted root mean square poloidal velocity 
(rms velocity) by 



I„,dj77(Vp-Vp) 

Lm dm 
Emdravn 



(19) 
(20) 



with the poloidal velocity vector Vp. 

The dependency of Vims on time is shown for all runs in Fig- 
ure[TT] The clouds set up in pressure equilibrium show first a shal- 
low rise and then stay roughly constant. From the density his- 
tograms (Figure [6}, the acceleration decreases when the densest 
cores are completely shredded. Hence, it seems that kinetic energy 
may be generated as long as dense cloud material may be squeezed 
radially and pushed out vertically. Since run R36 has a higher res- 
olution than run R34, the gas may be compressed further. More 
energy is therefore stored temporarily in the magnetic field, which 
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Figure 7. Close-up of the cloud core of run R63 at the beginning (top row) and at the end (bottom row) of the simulation. The logarithm of the density is shown 
on the left, the azimuthal magnetic field in the middle and the magnitude of the poloidal field with overlaid poloidal field vectors on the right. The coordinates 
are relative to the initial cloud centre. The slight upwards shift at the end of the simulation is due to an interaction with the boundary, which defines the end of 
the simulation. 
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Figure 8. Comparison of the density distribution of run R34 (left) to run R36 (double resolution, middle) and run R47 (double spatial scales, right) at the 
common snapshot time of 4.58 morb. 
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Figure 9. Comparison of the density distribution of run R60 (middle, smaller grid but otherwise identical to R63) to run R62 (half resolution, left) and run R61 
(50 per cent higher resolution, right) at the common snapshot time of 48.58 morb. 



ultimately results in higher velocities. Run R47 has twice the spa- 
tial scales but is otherwise identical to run R34. The net force per 
unit volume has not changed. Since the initial densities are also 
the same, the accelerations in each cell are the same, too. There- 
fore, the turbulent velocities increase in a very similar way up to 
about 35 morb. Since the velocities are the same, but the length 
scale is twice as long in run R47, there is still room for more com- 
pression in simulation R47 at this time. Consequently, the turbulent 
velocity increases further, until the resolution limit is reached near 
57 morb. One may understand the final turbulent velocity in terms 
of the total available space for compression: The gravitational and 
centrifugal forces act at similar strength in every cloud cell. The 
radiative force is concentrated on the illuminated surface. The net 
compression force is roughly F = (1 — V'^)Fq, where V is the ini- 
tial rotation velocity of the cloud in Kepler units. The total energy 
transfered to the cloud in the limit of complete linear compression 
would therefore be: 
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Figure 10. Density histograms for runs R60 (solid black) R61 (dotted yel- 
low, 50 per cent higher resolution) and R62 (dashed green, half resolution) 
at 48.58 morb. The high density part disappears at low resolution. 



where the numerical factor is due to the average cloud thickness of 
KRf.\^/2. If this energy would be completely used for accelerating 
the cloud, a turbulent velocity of: 
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would result for V = 0.85. As may be seen from Figure [TT] the 
velocities are somewhat higher than this estimate. This is likely 
due to the increased radiative force during the lateral expansion of 
the clouds. 

The clouds not set up in pressure equilibrium (R35 and R63) 
suffer initially a very quick increase of the rms velocity, as might 
have been expected. They evolve quickly towards pressure equilib- 
rium. Then, as in the initial pressure equilibrium runs, the turbulent 
velocities change only gradually. Several curves show an upturn to- 
wards the end, which is due to the beginning interaction with the 
grid boundary. 



6 DISCUSSION 

We have investigated whether magnetic fields might contribute to 
the stability of BLR clouds. Stability is an issue in models with- 
out magnetic fields (compare section[Tll. We first combine literature 
data to get an idea of the magnitude of magnetic and other forces 
that might be relevant in this context. We find that reverberation 
mapping data is consistent with the bound cloud picture for a rea- 
sonable column density. Yet, the forces of radiation and gravity are 
comparable. We argue that the magnetic energy density, which sets 
the magnitude of the magnetic pressure and tension forces, should 
be at least of the same order of magnitude as the radiation pressure, 
in order to be at all able to contribute to the stability of the cloud. 
Estimating these parameters from the literature turned out to be 
complicated by the fact that the flux of radiation inferred via rever- 
beration mapping studies leads to a rather high photoionisation pa- 
rameter. One could take this as indicati on that there are no c louds at 
all and cloudless models like the one o f lMurravetal] ( ll995b should 
be preferred. Yet, for cloud densities on the high side of the allowed 
values, which would also mean that the typical BLRS would more 
likely be gravitationally bound, or when adopting other changes 
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Mass weighted turbulent velocity over time 
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Figure 11. Mass-weighted root mean square velocities for runs R34 (black 
plus signs), R35 (blue stars), R36 (green squares), R47 (red diamonds) and 
R63 (orange triangles) over time. Only the poloidal velocity component is 
taken into account. 



of the relevant parameters within the range allowed by the uncer- 
tainties, the photoionisation parameter may be brought down to a 
reasonable value. 

A gravitationally bound cloud model entails the problem that 
the radiation leads to a surface force whereas gravity and centrifu- 
gal forces are body forces. Hence, an equilibrium configuration 
with these forces alone may not be found for the aforementioned 
assumptions. As all these forces act only radially away or towards 
the SMBH, an isotropic thermal pressure that was to stabilise the 
cloud against the compressive forces would necessarily also cause 
lateral expansion. Magnetic forces are anisotropic and are therefore 
in principal suited to counteract the compression and stabilise the 
clouds. We have performed MHD simulations with different mag- 
netic field configurations to test this idea. 

With an azimuthal magnetic field, only, the dynamics is dom- 
inated by the open tube of toothpaste mechanism: radial com- 
pression is followed by vertical expansion. Kelvin-Helmholtz and 
column density instabilities together with repeated, but smaller, 
tube of toothpaste squeezing events transform the initially homo- 
geneous cloud quickly into a filamentary system. This result is in 
good agreement with the unstable pancake picture (Mathews 1982, 
Il986h . In fact, in this case the effect of the magnetic field is just an 
additional pressure. Hence, the result remains essentially the same 
as in the pure hydrodynamic case. The cloud fragments are shred- 
ded down to the resolution level. Mixing with the pressure sup- 
ported ambient gas reduces the rotational support, and the clouds 
fall inwards. Additionally, the atmosphere is not entirely stable and 
features some inwards flow due to the cooling. The timescale for 
this process is given first by the cloud compression timescale: 



4R, 



1/2 



(23) 



4 months ( — 



-1/2 



V9200/xpc ) V3/tpc ) 



1/2 



where g is the gravitational acceleration, and the evaluation is again 
for an initial rotation velocity V — 0.85. Our numerical simula- 
tions show that clouds with an azimuthal magnetic field, only, 
are dispersed to the resolution level, and hence rapidly dissolve 
in the ambient medium within about tr, which is also similar to 



the Kelvin-Helmholtz timescale for the simulated cases. We show 
that the dispersal is even faster for higher resolution. The mix- 
ing, on the other hand, depends on the ambient density, via the 
Kelvin-Helmholtz timescale. It is therefore expected to be slower 
for lower ambient densities. Not much is known about the den- 
sity in the inter-cloud gas. Assuming that the density would be 
lower by a factor of lO'*, would increase the Kelvin Helmholtz 
timescale by a factor of 100, to the range of the orbital timescale. 
This would lead to ext remely spread ou t pancake clouds, similar to 
the ones envisag edbv lMathewsl ( ll986l) . We have also tried a fac- 
tor of ten lower inter-cloud density than reported here. In this case, 
the cloud spreads quite quickly over the grid boundary probably 
due to the diminished resistance of the ambient gas, whereas the 
density histograms showed more dense gas at comparable times, as 
expected. In three dimensions, such pancake clouds would also face 
the ram pressure, which might also lead to compressed cloud heads 
with cometary tails, and a stability timescale of order the orbital 
timescale (Math ewsll 19861) . 

The runs with azimuthal magnetic field, only, underline the 
stability problem, and demonstrate clearly that our simulations are 
able to capture the destruction process adequately. Our runs with a 
helical magnetic field result in a much more stable cloud configura- 
tion, which survives the destruction process for the simulation time. 
The cloud corresponds to a helical filament in azimuthal direction. 
The core of this filament is dominated by the azimuthal field, pro- 
ducing an overpressure of about 30 per cent. This is balanced by the 
magnetic tension force of the poloidal field surrounding the core. 
The stability of the configuration results from the nature of the ten- 
sion force which is inversely proportional to the curvature radius 
of the field lines. A radial compression straightens the field lines 
vertically (less radial tension force) and bends them more sharply 
at the vertical ends (increased vertical restoring tension force). A 
similar restoring force is obtained for compressions in every other 
direction. This configuration is very similar to the magnetic field 
used in tokamaks. One reason why it is employed in this context is 
precisely its MHD stability. 

In order to isolate the effects of the magnetic field, we have 
used a negligible thermal cloud pressure. In our simulations we 
have realised this via a minimum temperature of 1000 K. The sim- 
ulations would not have been any different had we used 10"* or even 
10^ K, as the magnetic pressure would still be by far dominant in 
the clouds. We had initially decided for a low minimum temper- 
ature to check if the induced dynamics would lead to any signifi- 
cant increase of the cloud temperature via shock heating. The bal- 
ance between shock heating and radiative cooling was found by 
iKrause & Alexanded ( |2007|) in their simulations of clouds in jet 
cocoons to produce a peak in the temperature histograms around 
10^ K. We do not find this effect here: Even for the late phases of 
R34, where the cloud is already essentially dispersed, the by far 
major part of the filaments remain close to the minimum tempera- 
ture. Therefore, we do not find a heating mechanism which could 
compete with photoionisation. 

Cloud stability in the limit of negligible magnetic fields has 
been studied extensively in the past (compare section [TJ. We have 
studied here the case of a dominant magnetic field compared to 
the thermal energy density, and separated the effects of magnetic 
pressure and magnetic tension. The azimuthal field behaves like a 
pressure in our axisymmetric simulations. Hence, we expect that 
we could in principle easily accommodate for a thermal pressure 
comparable to the magnetic one by replacing some or all of the 
magnetic pressure due to the azimuthal field by a corresponding 
thermal pressure. Details would of course depend on the effective 
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equation of state, as the thermal pressure is affected by heating and 
coohng processes. 

Another question one has to ask is how likely it would be to 
find such a helical configuration around real AGN. If the inter-cloud 
density would be comparable to the one assumed in our simula- 
tions, all unstable configurations would be shredded and dissolved 
within months. It is then quite difficult to replace the clouds quickly 
enough. If they would for example be launched from the accretion 
disc by a wind with velocities comparable to ones observed in the 
BLR, it would take a similar amount of time to be lifted to a signifi- 
cant altitude abov e the disc. Yet, the BLR cannot be a thin disc (e.g. 
IOsterbrock|[l978h . Our helical field simulation has started quite far 
away from the final equilibrium, yet it has reached a stable con- 
figuration. It seems therefore possible that many different kinds of 
initial conditions eventually find an equilibrium solution. It is be- 
yond the scope of this study how many stable filaments may be 

produced in given conditions of the interstellar medium. 

Regarding the absolute magnitude of the magnetic field. iRee^ 

( [l987i) has already pointed out that magnetic field values of a few 
Gauss, as required here, are expected in BLRs due to relativistic 
winds, accretion flows or accretion driven winds. He also mentions 
the possibility of magnetic fields within the clouds if they are cre- 
ated by the thermal instability. We note that a similar value for the 
magnetic field has recently been estim ated from their possib le ef- 
fects on the polarisation of the Ha line jSilant' ev et al.ll20I2h . 

What level of internal turbulence do our simulations predict? 
We have argued that the available energy for turbulent motions is 
limited by the cloud compression and the Alfven speed. Since the 
initial cloud volumes are unknown, we cannot make a quantitative 
prediction from this argument. If the BLR cloud population would 
be dominated by unstable clouds, being constantly shredded and 
replaced, we would predict turbulent velocities of about a third of 
the magnetosonic speed. If the radiation pressure would indeed be 
matched by the magnetic pressure in the cloud, the magnetosonic 
speed would be of order 1000 km/s, and the predicted internal ve- 
locities consequently a few hundred. For stable clouds of the heli- 
cal type as shown above, we would expect even less in the cloud 
co res. Hence, the turbu l ent ve locities of about 1000 km/s derived 
bylKollatschny & Zetzl|j 201lh cannot be reached in this way. 

The cometary clouds observed bv lMaiolino et al.l ( l2010h . with 
dense cores and filamentary tails are compatible with our simu- 
lation results. Taking the orbital motion into account, the filamen- 
tary tails would be dominantly elongated in the azimuthal direction, 
as probably also in the observations. If these cloud s would indeed 
be des troyed within a few months, as estimated bv iMaioIino et al.l 
l l2010h . this would fit exactly with the destruction timescale of the 
radiative destruction mechanism discussed here. 

In the literature one may find also BL R models that involv e 
no clouds at all, e.g. the disc wind model of iMurrav et alj jl995l) . 
Here, the BLR forms the base of a wind accelerated by radiative 
and thermal pressure. The advantages and disadvantages of models 
with a nd without clouds have been nicely summarised bv iNetzeJ 
( l2008h : Cloud models provide a significantly better match to the 
emission line structures, while cloudless models obviously avoid 
the confinement problem altogether. 

We have addressed the meridional stability problem, only. 
While this is certainly one of the major issues in the BLR cloud 
stability problem, we have not shown that the cloud would also 
survive 3D effects. From the MHD point of view, kink instabilities 
are likely to occur, which would cause bends in the cloud filament 
along its major axis. Further, the clouds are unlikely to be extended 
along the entire azimuthal angle as necessarily implied in our axi- 



symmetric simulations. In reality, the clouds could still be elon- 
gated substantially in the azimuthal direction, if the ends would 
be anchored in a thinner disc structure, similar to coronal loops 
on the surface of the sun. In simple magnetic field configurations 
with closed field lines within the cloud, the field could of course 
be transformed via reconnection, which would likely lead to cloud 
splitting as each closed field region could support a cloud on its 
own. An azimuthally extended cloud with the azimuthal field clos- 
ing at the ends and returning through some part of the cloud would 
be expected to expand in the polar direction and form something 
like a closed ring. For any non-axisymmetric clouds, also the ram 
pressure would be significant. For our cloud setups the azimuthal 
ram pressure would be comparable to the magnetic pressure in the 
cloud. If one can find magnetic field configurations that would also 
stabilise against this effect is beyond the scope of the present inves- 
tigation. 

In order to reach pressure equilibrium with the clouds, also 
the inter-cloud medium should be magnetised on a similar level. 
This might be accessible to Faraday rotation. The expected rotation 
measure would be of order 

RM = 8. 12x10** /^;1e6B•d/radm"^ (24) 
JO 

where the path-length d/ is measured in units of 10^^ pc, B in Gauss 
and the electron density in units of 10^ cm^^ . As mentioned above, 
the density is uncertain, and could also be a few orders of magni- 
tude lower. For a turbulent medium, one would therefore predict 
that signals below ~ lO'^ Hz {X > 30/im) would be depolarised 
jtCrause. Alexander. & Boltonll2007h . This prediction might there- 
fore be tested in the future by infrared polarimetry, if one is sure to 
observe emission from a jet base. 

For our simulations, a change in the geometry of the magnetic 
field in the inter-cloud region should not matter much, because the 
thermal energy density dominates there over the magnetic one. As 
long as the forces are set up close to equilibrium, we would expect 
an outcome in between runs R34 (initial pressure equilibrium) and 
R35 (no initial pressure equilibrium). In reality, the situation might 
of course be more complex, and the inter-cloud medium might 
be magnetically dominated and at the same time have a non-zero 
poloidal component. It would then be important if the magnetic 
field would be topologically connected to the clouds. In that case, 
the cloud gas could escape along the field lines. It might also hap- 
pen that regions of space with a favourable magnetic field config- 
uration might protect new clouds which are just forming from the 
thermal instability. Yet, these complexities are beyond the scope of 
the present discussion. 

Summarising, we find some indication that bound BLR clouds 
might indeed be stabilised by the magnetic field against collapse 
due to opposing radiative and gravitational forces. In this case, the 
magnetic field would have to have a poloidal component. 



7 CONCLUSIONS 

Gravitationally bound clouds facing strong radiation pressure are 
unstable in the purely hydrodynamic case because radiative and 
gravitational forces compress the clouds radially, whereas the ther- 
mal pressure acts isotropic. They may be more stable if they are 
significantly magnetised. In particular, we find in axisymmetric 
magnetohydrodynamic simulations with a prescription for the ra- 
diation pressure that the magnetic tension force produces more sta- 
ble clouds, while in a situation where the geometry of the magnetic 
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field is such that its effect is only analogous to an additional pres- 
sure, the clouds are similarly unstable as in the hydrodynamic case. 
In order to be effective, in the BLRs of a discussed sample of AGN, 
the magnetic field strength should be of order a few Gauss, accu- 
rate to about an order of magnitude and constrained by the con- 
dition that magnetic, radiative and gravitational forces should be 
comparable. 
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